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Abstract. We consider the heavy quark limit of Coulomb gauge QCD, with the truncation of the Yang-Mills sector to include 
only (dressed) two-point functions. We find that the rainbow-ladder approximation to the gap and Bethe-Salpeter equations is 
nonperturbatively exact and moreover, we provide a direct connection between the temporal gluon propagator and the quark 
confinement potential. Further, we show that only bound states of color singlet quark-antiquark (meson) and quark-quark 
(SU(2) baryon) pairs are physically allowed. 
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Coulomb gauge is an ideal choice for investigating the 
confinement phenomenon. In this gauge, there is an ap- 
pealing picture of confinement: the Gribov-Zwanziger 
scenario |U], whereby the temporal component of the 
gluon propagator provides a long-range confining force 
whereas the transverse spatial components are infrared 
suppressed (and therefore do not appear as asymptotic 
states). Confinement implies the existence of a nonper- 
turbative scale (the string tension), and this naturally 
leads to the question: is there a simple connection be- 
tween the (gauge fixed) Green's functions of Yang-Mills 
theory and this physical scale? To answer this ques- 
tion, we propose to study here the full (nonperturbative) 
QCD Bethe-Salpeter equation in Coulomb gauge within 
a heavy quark mass expansion at leading order. In this 
truncated system, we will use the results inspired by the 
lattice for the explicit Green's functions of the Yang- 
Mills sector and in addition we will employ the Slavnov- 
Taylor identity for the quark-gluon vertex. 

We begin by considering the explicit quark contribu- 
tion to the full QCD generating functional (unless oth- 
erwise specified, the Dirac and color indices in the fun- 
damental representation are implicit and we follow the 
conventions from (2D): 
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Sym represents the Yang-Mills part of the action and 
the temporal and spatial components of the covariant 
derivative (in the fundamental color representation) are 



given by 

D = d -igT a a a (x), D = V + igT a A a (x), (2) 

where A and a refer to the spatial and temporal compo- 
nents of the gluon field, respectively. 

The full quark field q (q is the conjugate field, and %, % 
are the corresponding sources) is decomposed according 
to the heavy quark transformation: 

q{x)=e-" nx °[h(x)+H(x)} 7 
1 + 7° 



h(x 



mi.x ; 



q(x), H(x)=e'"^l-^q(x)(3) 



(similarly for the antiquark field). We now insert this de- 
composition into the generating functional Eq. £[), in- 
tegrate out the //-fields, and make an expansion in the 
heavy quark mass, similar to Heavy Quark Effective The- 
ory [HQET] ( 3]. At leading order, we get the following 
expression: 

Z\X,X} = l^expj* Jd 4 xh(x) [id 0x + gT a O a ]h{x) 

+ ijd 4 x [e- l ^x{x)h{x) + e im *h{x)x{x)] + i^ym 

+0(l/m). (4) 

In the above, we have retained the full quark and anti- 
quark sources (unlike HQET). This means that we can 
use the full gap and Bethe-Salpeter equations of QCD 
replacing, however, the kernels, propagators and vertices 
by their leading order expression in the mass expansion. 
Also, notice that the spin degrees of freedom have de- 
coupled from the system. 



In full Coulomb gauge QCD, the quark gap equation 
is given by [J 
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(r's denote the various proper functions, W denotes 
propagators, see [4]). The gap equation is supplemented 
by the Coulomb gauge Slavnov-Taylor identity 01: 
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where k\ + k 2 + k 3 = 0, qo is an arbitrary energy injec- 
tion scale (arising from the noncovariance of Coulomb 
gauge), Tec is the ghost proper two-point function, V-qpcq 
and Tqxcq are ghost-quark kernels associated with the 
Gauss-BRST transform (see also 

Now, as a consequence of the (Coulomb gauge) de- 
composition, Eq. (O and under the assumption that 
the pure Yang-Mills vertices may be neglected (re- 
taining, however, the dressed gluon propagator), the 
Dyson-Schwinger equation for the nonperturbative spa- 
tial quark-gluon vertex furnishes the result that TqqA ~ 
0(1 jm) (see 10] for a complete discussion and justi- 
fication of this truncation). Similarly, the ghost-quark 
kernels can be neglected. Thus, under our truncation 
scheme, the Slavnov-Taylor identity reduces to 
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This is then inserted into Eq. (O, together with the tree- 
level quark proper two-point function 



and the tree level quark gluon vertex 



(8) 
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that follow from the generating functional Eq. (|4j. The 
general form of the nonperturbative temporal gluon prop- 
agator is given by @1: 
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1 This expression is in the second order formalism and can be directly 
inferred from the result obtained within the first order formalism 141. 



Lattice results [7] motivate that the dressing function 
Daa is largely independent of energy and likely to be- 
have as 1/P for vanishing k 2 (the explicit form of D aa 
will only be needed in the last step of the calculation). 
Putting all this together, we find the following solution 
to Eq. (O for the heavy quark propagator: 



Wqq(k ) = 
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where the (implicitly regularized, denoted by "r") con- 
stant is given by 
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with the Casimir factor C F = (N 2 - l)/2N c . When solv- 
ing Eq. (0, the ordering of the integration is set such that 
the temporal integral is performed first, under the con- 
dition that the spatial integral is regularized and finite. 
Inserting the solution Eq. (fTTT i into the Slavnov-Taylor 
identity, we find that the temporal quark-gluon vertex re- 
mains nonperturbatively bare: 
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Note that the propagator Eq. (fTTT i has a single pole in 
the complex fcrj-plane (due to the mass expansion) and 
therefore it is necessary to explicitly define the Feynman 
prescription. From Eq. (fTTT i it then follows that the closed 
quark loops (virtual quark-antiquark pairs) vanish due to 
the energy integration, which implies that the theory is 
quenched in the heavy mass limit: 
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We also emphasize that the position of the pole has no 
physical meaning since the quark can never be on-shell. 
The poles in the quark propagator are situated at infinity 
(thanks to J? r as the regularization is removed) meaning 
that either one requires infinite energy to create a quark 
from the vacuum or, if a hadronic system is considered, 
only the relative energy (derived from the Bethe-Salpeter 
equation) is important. 

The solution for the antiquark propagator reads: 
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Notice the assignment of the Feynman prescription, sim- 
ilar to Eq. ( fTTl ) - see also the discussion in yfl. This has 
the important consequence that the Bethe-Salpeter equa- 
tion for the quark-antiquark states has a physical inter- 
pretation of a bound state equation (see below). The cor- 
responding vertex reads: 
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Let us now consider the full homogeneous Bethe- 
Salpeter equation for quark-antiquark bound states: 

r(p;P) a /3 = ~(2^j4 / dkK aP;S r (P' k ' P ) 

x [W m {k + )T{k;P)W lq {k-)] rS (17) 

where k± (similarly for p±) are the momenta of the 
quarks (with the notations from [2]), P is the 4- 
momentum of the bound state (assuming that a solu- 
tion exists), K represents the Bethe-Salpeter kernel and 
r is the Bethe-Salpeter vertex function for the particu- 
lar bound state under consideration. Further, we explic- 
itly identify the antiquark contribution, i.e. Wq q (k-) — 
—Wqq{—k-) (also in the kernel). 

When constructing the Bethe-Salpeter kernel, we 
use the fact that the temporal integration performed 
over multiple quark propagators with the same relative 
sign for the Feynman prescription vanishes (similar to 
Eq. (O, but in this case the terms originate from inter- 
nal quark or antiquark propagators) and hence, the kernel 
reduces to the ladder truncation Q2[| : 

K a p; Sr (p,k) = r qqaar w^(k)r^ ap5 . (is) 

We now insert the nonperturbative results for the prop- 
agators and vertices, Eqs. (II 11131 11511 6b and take the 
form, Eq. dlOl . for the temporal gluon propagator. Fur- 
ther, we notice that that the Bethe-Salpeter equation is in- 
dependent of the relative quark energy and hence we can 
perform the temporal integration over the quark propaga- 
tors. The energy integration over the quark and antiquark 
propagators now leads to (unlike Eq. (TBI): 



J-oo [k° + - m - J r + is] [k°_ - m + J? r - ie] J> r 

By using the expression Eq. ([T2"i i for J^ r , and Fourier 
transforming to coordinate space, we find the following 
simple energy solution for the pole condition of the 
Bethe-Salpeter equation Eq. ( fTTI i: 

(20) 

where Cm arises from the color structure and is yet to be 
identified (r is not assumed to be a color singlet): 

[T a mT a } ap =C M r ap (y). (21) 

Because the total color charge of the system is con- 
served and vanishing [8], a single quark (or antiquark) 
cannot be prepared in isolation. Thus, the bound state 
energy Pq can only be either confining for large sepa- 
rations, i.e. increase linearly with the separation between 
the quark and antiquark, or be infinite when the hypothet- 
ical regularization is removed (so that the system cannot 



be physically created). If the temporal gluon propagator 
dressing function is more infrared divergent than l/|d)|, 
then 

Cf = Cm (22) 

is required such that the spatial integral is convergent. 
This gives the condition 

r ay(y) = <S ar r(y), (23) 

which means that the quark-antiquark Bethe-Salpeter 
equation can only have a finite solution for color singlet 
states and otherwise the energy of the system is diver- 
gent. Assuming that in the infrared D aa =X/a> 2 (as in- 
dicated by the lattice) where X is some combination of 
constants, then from Eq. (120b with the condition Eq. (l22i i 
we find 

P = a\y\ = ^^\y\ + ^(l/m). (24) 

The above result is that there exists a direct connection 
between the string tension a and the nonperturbative 
Yang-Mills sector of QCD at least under the truncation 
scheme considered here. 

A similar calculation performed for the diquark Bethe- 
Salpeter equation shows that the diquarks are confined 
for N c = 2 colors, corresponding to the SI/ (2) baryon, 
and otherwise there are no (finite) physical states. 
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